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V i Properties and mechanism of the magnetic phase transition of the perovskite-type Ti 

H ■ 

oxides, which is driven by the Ti-O-Ti bond angle distortion, are studied theoretically 



by using the effective spin and pseudo-spin Hamiltonian with strong Coulomb repulsion. 
It is shown that the A-type antiferromagnetic(AFM(A)) to ferromagnetic(FM) phase 
*^ . transition occurs as the Ti-O-Ti bond angle is decreased. Through this phase transition, 

the orbital state is hardly changed so that the spin-exchange coupling along the c-axis 
changes nearly continuously from positive to negative and takes approximately zero 
at the phase boundary. The resultant strong two-dimensionality in the spin coupling 
. causes a rapid suppression of the critical temperature as is observed experimentally. 
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§1. Introduction 

Perovskite-type Ti oxide RTiQ^R being a trivalent rare-earth ion) is a typical Mott- 
Hubbard insulator^. Ti 3+ has a 3d 1 configuration, and one of the three-fold ti g orbitals is 
occupied at each transition-metal site. The crystal structure is an orthorhombically distorted 
perovskite (GdFeOs-type distortion) whose unit cell contains four octahedra as shown in 
Fig. 1. The magnitude of the distortion depends on the ionic radii of the R ion. With a 
small ionic radius of R ion, the lattice structure is more distorted and the bond angle is 
decreased more largely from 180°. The bond angle can be controlled by use of the solid- 
solution systems Lai-yY^TiOa or in RTi0 3 , varying the R ions. Especially, with varying the 
Y concentration in Lai^Y^TiOs, we can control the bond angle almost continuously from 
156° (y = 0) to 140° (y = 1). In YTiOs, a <i-type JT distortion has been observed where 
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Fig. 1. GdFeC>3-type distortion. 



the longer and shorter Ti-0 bond lengths are ~2.08A and ~2.02A, respectively® . Although 
the difference between the longer and shorter bond length is relatively small, LaTi03 also 
shows a G?-type JT distortion. Recently, electronic and magnetic phase diagrams have been 
investigated intensively as functions of the magnitude of a Ti-O-Ti bond angle distortion 
00 Si), in the less distorted or a La-rich(y < 0.6) region, the system shows AFM ground 
state. In particular, LaTi03 (y = 0.0) shows a G-type AFM(AFM(G)) ground state with a 
magnetic moment 0. With increasing the Y-concentration or varying the R site with 

smaller size ion (decrease of the Ti-O-Ti bond angle), the Neel temperature (Tn) decreases 
rapidly and is suppressed to almost zero, and subsequently a FM ordering appears. This 
rapid decrease of Tn is hardly explained by the conventional models. The origin is one of 
the issues of interest. In the relatively distorted or Y-rich region, the system shows a FM 
ground state. In YTiOs(?/ = 1.0), the value of the magnetic moment is 0.84/i B i ) . This 
ferromagnetism is hardly explained by a simple single-band Hubbard model and requires 
to consider the ci-level degeneracy. Recent model Hartree-Fock studies have succeeded in 
reproducing the magnetic structures of the both end compoundsi'0). The nature of the 
magnetic phase transition is, however, not sufficiently clarified. Besides, the spin and orbital 
states realized in the moderately distorted region are issues of interest. In this paper, we 
study the properties and mechanism of the phase transition by focusing on the region near 
the phase boundary. 

§2. Formalism 

We start with the multi-band d-p model in which the full degeneracies of Ti 3d and O 
2p orbitals as well as on-site Coulomb and exchange interactions are taken into account. 
In this Hamiltonian, the effects of the GdFeO"3-type distortion are considered through the 
d-p transfer integrals which is defined by using the Slater-Koster's parameters Vp dn ,Vp da ,Vp piT 
and V^, po 0. The effects of the <i-type Jahn-Teller(JT) distortion are also considered. The 
magnitude of the distortion is expressed by the ratio [V pda ./V pdo ] 1//3 , here V pda and V pda are the 
transfer integrals for the shorter and longer Ti-0 bonds. In order to reveal the origin of the 
rapid suppression of Tn, we focus on the situation near the phase boundary between AFM 
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and FM phases. The value of the ratio [V^/V^.] 1 / 3 is fixed at 1.030, which is expected to be 
realized near the phase boundary under the assumption of the linear increase as a function 
of the bond angle from 1.036(YTiO 3 ) to 1.00(LaTiO 3 ). Under the JT distortion, the t 2g 
level-splitting energy A t . 2g is estimated as 0.050eV by using Slater-Koster's relations. Since 
the t 2g level-splitting due to the spin-orbit interaction is sufficiently small in comparison with 
Af 2s , we neglect the interaction through the present calculations. By integrating over the O 
2p orbital degrees of freedom, we derive the effective Hubbard model which includes only 
Ti 3d orbital degrees of freedom. 

Under the cubic-type crystal field, a five-fold degeneracy of 3d orbitals is lifted to two-fold 
higher levels e g and three-fold lower levels t 2g . Moreover, in the d-type JT distortion, the 
degeneracies of the e g and t 2g orbitals are lifted. Although the e g level is uniformly lifted to 
a higher level 3z 2 — r 2 and a lower level x 2 — y 2 at each site, the ways of t 29 Tevel splitting 
are different between sites 1, 3 and sites 2, 4 (see Fig. 1). In sites 1 and 3, the y-axis is 
elongated and consequently, the t 2g level is lifted to lower xy and yz levels and a higher 
zx level. On the other hand, in sites 2 and 4, the x-axis is elongated so that it is lifted to 
lower xy and zx levels and a higher yz level. Let us represent the five 3d orbitals xy, yz, 
zx, x 2 — y 2 and 3z 2 — r 2 by energy-level indices 1, 3, 2, 4 and 5, respectively at sites 1 and 
3, and by indices 1, 2, 3, 4 and 5 at sites 2 and 4. As a result, in the insulating d 1 systems 
under a JT distortion, one of the two-fold degenerate lower t 2g orbitals is occupied at each 
site. Especially, in the case of the d-type JT distortion, either xy or yz orbital is occupied 
at site 1 and 3, xy or zx at site 2 and 4. 

Based on the above discussion, the multi-band Hubbard Hamiltonian derived from the 
multi-band d-p model has a form, 



H = tl d + H tM + -Hon-site l-U 

with 

H-d = £di,md\^ m ^di trrlt<T , (2) 
i,m,cr 

^tdd = tim,i'm'd'i,m,tjd'i' ,m' ,<j + h.C, (3) 
i,m,i' ,m' ,cr 

-^on-site = H u + H u i + Hj + Hji, (4) 

where d\ m u is a creation operator of an electron with spin a(= \, |) in a 3d level m at Ti 
site i. By H™ 11 , we express the level energies of Ti 3d orbitals under the influence of the 
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crystal fields with 

s dl for m = l,3, 



£di + At 2s for m = 2, 
£<« + A z 2_ y2 for m = 4, 
e<8 + A 3z 2_ r 2 for m = 5. 
Here, m = 1, 3 are lower t 2 g levels, m = 2 is a higher t 2 g level and m = 4 and m = 5 are lower 
and higher e 9 levels, respectively. The \ 2g , A x ' 2 -y 2 and A 3z 2_ r 2 denote the level-splitting 
energies measured from lower t2 g level. It should be noted that the same indices of energy 
levels at different sites do not necessarily correspond to the orbitals with the same symmetry. 
H^ is a d-d super-transfer term and H on _ site represents on-site d-d Coulomb interactions. 
The -f/on-site term consists of the following four contributions, 

H u = ^ud\ m] d itmA d\ mi d i)m)[ , (6) 

i,m 

H u ' = U d>i, m ,adi,m,ad iim /ydi jm iy, (7) 

i,m>m' ,<j,<t' 

Hj = jd imo .di trn i^ <J d i , m , IJ ,di tm jji, (8) 

i,m>m' a,cr' 

Hj' = f^i,m,^i,m',ldi, m ,[di,m',l, (9) 

where H u and H u > are the intra- and inter- orbital Coulomb interactions and Hj and Hf 
denote the exchange interactions. The term Hj is the origin of the Hund's rule coupling 
which strongly favors the spin alignment in the same direction on the same atoms. These 
interactions are expressed by using Kanamori parameters, u, u', j and j' which satisfy the 
following relations^©, u = U + yj, v! = u — 2j and j = j'. Here, U gives a 
magnitude of the multiplet-averaged d-d Coulomb interaction. The charge-transfer energy 
A, which describes the energy difference between occupied O 2p and unoccupied Ti 3c? levels, 
is defined by using U and energies of the bare Ti 3d and O 2p orbitals e° d and e p as follows, 

A = e° d + U - e p , (10) 

since the characteristic unoccupied 3d level energy on the singly occupied Ti site is e° d + U. 
The values of A, U and V p da are estimated by photoemission spectraB @>. We take the values 
of these parameters as A = 6.0eV, U = 4.0eV, V v da = — 2.0eV and j = 0.74eV throughout 
the present calculation. The ratio V^/V^- is fixed at —2.17 and V ppa and V PJm at 0.60eV 
and —0.15 eV, respectively©. 

Starting with the multi-band Hubbard Hamiltonian, we can derive an effective Hamilto- 
nian in the low-energy region on the subspace of states with singly occupied t2 9 orbitals at 
each transition-metal site by utilizing a second-order perturbation theory. The states of 3d 
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Fig. 2. Relative energy as a function of the Ti-O-Ti bond angle. 



electron localized at the transition-metal sites can be represented by two quantum numbers, 
the z-component of the spin S z and the number of the occupied orbitals. When one of the 
two-fold lower t 2g orbitals is occupied at each site, we can describe the electronic states using 
a spin-1/2 operators, which we call the pseudo-spin r. We follow an approach similar to the 
well-known Kugel-Khomskii formulation§>@0§). We express the 3d electron operators in 
terms of S and r to arrive at the effective spin and pseudo-spin Hamiltonian; 

H cS = H™ u + H t2g + H eg . (11) 

The first term H™ R is obtained from the zeroth-order perturbational processes. The second 
term H t2g is obtained from the second-order perturbational processes whose intermediate 
states contain only ^-orbital degrees of freedom. The third term H eg is obtained from the 
second-order perturbational processes whose intermediate states contain e 9 -orbital degrees 
of freedom. 

§3. Results 

We have calculated the total energies of various spin and orbital configurations by applying 
a mean-field approximation to the effective spin and pseudo-spin Hamiltonian. Near the 
AFM-FM phase boundary, the AFM(A) and FM solutions with a certain kind of orbital 
orderings are stabilized. Figure 2 shows that the AFM(A) to FM phase transition arises by 
decreasing the Ti-O-Ti bond angle (increasing the GdFe03-type distortion). In the AFM(A) 
and FM phases, yz, zx, yz and zx orbitals are predominantly occupied among the two-fold 
degenerate lower t 2g orbitals at site 1, 2, 3 and 4, respectively ((yz, zx, yz, zx)-tjpe orbital 
order). We can specify the orbital states realized in the AFM(A) and FM solutions by using 
the angle #afm(A) and 9^m as, 

sitel; cos9 x \xy > + sm9 x \yz >, 
site2; cos 9 x \xy >+ sin 9 x \zx >, 



5 



site3; — cos6 x \xy > + sin 6 x \yz >, 

site4; — cos 9 X \ xy > + sin 6 x \zx >, (12) 

where x = AFM(A), FM. In Fig. 3, the angles for the AFM(A) and FM solutions are plotted. 
The difference between the #afm(A) and 9 FM is very small, especially in the more distorted 
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Fig. 3. Orbital states of the FM- and AFM(A) solutions. 



region or near the phase boundary. This means that the way of the orbital ordering hardly 
changes through the magnetic phase transition. Then the AFM(A) to FM phase transition 
is identified as the transition where the sign of the spin exchange interaction along the c-axis 
is changed from negative to positive while that in the a6-plane is constantly negative. The 
constant FM coupling in the afr-plane under the (yz, zx, yz, zx)-type orbital state can be 
easily understood. In the afe-plane, the neighboring orbitals (yz and zx) are approximately 
orthogonal to each other. Hence, the FM spin configuration is favored through Hund's 
coupling interaction. However, the emergence of the FM phase is still controversial since the 
neighboring orbitals along the c-axis are not orthogonal but have the same symmetries. By 
considering the transfers from yz to neighboring 3z 2 — r 2 , this is understood schematically 
as follows. 

At this stage, we assume that the (yz, zx, yz, ,2x)-type orbital order is hardly changed 
between two phases. Let us consider the energy gain of an electron in the yz orbital at site 
1, which is caused by the second-order perturbational processes with respect to the transfers 
along the c-axis (i.e., the transfers between site 1 and site 3). In the large GdFeC>3-type 
distortion, the yz orbital at site 1 mainly hybridize with the yz and 3z 2 — r 2 orbitals at site 3 
along the z direction relative to the other orbitals. When the yz orbital at site 3 is occupied 
by an electron, the second-order perturbational energy gain of an electron in the yz orbital 
at site 1 depends on the spin configuration between site 1 and site 3. When the spins of 
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Fig. 4. Characteristic second order perturbational processes. 



electrons on site 1 and site 3 are antiparallel, the absolute value of the energy gain can be 
written approximately as follows(see Fig. 4 a)), 

f 2 f 2 

- + ^~^ • (13) 

Here, t\ represents the transfer between yz at site 1 and yz at site 3 and t2 represents that 
between yz at site 1 and 3z 2 — r 2 at site 3 and A 3z 2_ r 2 denotes the level-energy difference 
between 3z 2 — r 2 and yz orbitals. On the other hand, when the spins are parallel, transfer 
to the yz orbital is forbidden by Pauli's principle but the energies of the intermediate states 
in which two electrons occupy the different orbitals are reduced by the intra-site exchange 
interaction j( see Fig- 4 b)). Consequently, the absolute value of the energy gain can be 
written as 

h h + M . (14) 

U' + A 3z 2_ r 2 — j U' + A 3z 2_ r 2 (u' + A 32 2_ r 2) 2 

Therefore, the spin configuration between site 1 and site 3 is determined by the competition 

of following two energies, — and , , , ? 2 In Fig. 5, the values of these energies are 

plotted as functions of the bond angle. As the GdFe0 3 -type distortion increases, the indirect 
hybridizations between neighboring i 2fl orbitals are decreased and those between neighboring 
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Fig. 5. Characteristic energy gains due to the second order perturbation. 



t 2g orbitals and e g orbitals are increased. Consequently, the value of t\/u is decreased and 
that of t\j / (u' + A 32 2_ r 2) 2 is increased, resulting in crossing of two energies as the bond angle 
is decreased. Moreover, since the hybridization between t2 9 and O 2p orbitals becomes to have 
a cr-bonding character, the amplitudes of the t% and t\j / (u'+A 3z 2_ r 2) 2 are critically increased 
by the GdFe0 3 -type distortion. Based on the above discussions, we can well describe this 
system by the following Heisenberg model as far as the (yz, zx, yz, zx)-type orbital order is 
strongly stabilized and hardly affected by the change of the spin configuration, 

c a,b 
-f^Heis — <^Heis ^ ' + ^Heis 5Z &i ' (^) 

(m) (hi) 

with 

< 0, (16) 
J Hcis = 4(2 - — - -f? V2 ). (17) 

Here, Yl%j) denotes the summation over the neighboring spin couplings along the c-axis 
and Y^tj) ^ n ^ ne a ^-pl a ne. The AFM(A) to FM phase transition occurs by the change 
in the sign of </H eis . Moreover, within this model, since the value of J^ eis decreases from 
a positive value to a negative one continuously as the bond angle is decreased and takes 
zero at the phase boundary. The two-dimensional spin coupling is realized at the phase 
boundary. Consequently, the critical temperature at the phase boundary suppressed to zero 
according to Mermin and Wagner's theorem©. Therefore, we can conclude that a strong 
two-dimensionality in spin couplings is realized near the phase boundary. 

In summary, the possible scenario of the rapid decrease of Tn is as follows. The GdFe0 3 - 
type distortion increases the indirect hybridizations between neighboring t2 9 and e g orbitals. 
As a result, the occupancy of the orbitals directed along c-axis is strongly favored by energy 



gain due to the second-order perturbational processes with respect to the transfers in the 
^-direction. Hence, the (yz, zx, yz, zx)-tjpe orbital order is strongly stabilized almost inde- 
pendently of the spin configuration. Since the orbital structure is hardly changed through 
the magnetic phase transition, the spin-exchange interaction along the c-axis is decreased 
from positive value to negative one nearly continuously with decreasing the bond angle, 
and consequently becomes almost zero at the phase boundary while that in the a6-plane 
remains constantly ferromagnetic. This strong two-dimensionality at the phase boundary 
suppresses T N and T c critically. However, a slight change of the orbital state exist at the 
phase boundary, and it causes the first-order phase transition in this system. 
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